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Abstract. A new method of extracting diffractive parton distributions is presented which 
avoids the use of Regge theory ansatz and is in much closer relation with the factorization 
theorem for diffarctive hard processes. 



1. Introduction 

Diffractive parton distributions functions (DPDF) are essential ingredients in the understanding 
and description of hard diffractive processes. The factorization theorem for diffractive process 
in Deep Inelastic Scattering (DIS) [TJ [2] enables one to factorize the diffractive DIS cross- 
sections into the long-distance contribution parametrised by DPDF's, from the short-distance, 
pertubatively calculable, one. Although DPDF's encode non-perturbative effects of QCD 
dynamics, their variation with respect to the factorization scale is predicted by pQCD O 2]. 
Moreover the short distance cross-sections is the same as inclusive DIS [I] so that higher order 
corrections can be systematically evaluated. Due to the factorization theorem, DPDF's are 
predicted to be universal distributions in the context of diffractive DIS. Next-to- leading order 
predictions based on DPDF's have been found to well describe DIS diffractive dijet cross- 
sections 0(6], thus confirming factorization. On the contrary, the issue of factorization is still 
debated in diffractive photoproduction of dijets since the HI [7] and ZEUS [8j collaborations do 
report conflicting results. 

DPDF's give the joint probability for a parton in a proton of four-momentum P to initiate a hard 
scattering keeping the proton intact with a four-momentum P' in the final state. Diffractive 
parton distributions are a special case of extended fracture functions [1] in the kinematic range 
xjp < 10~ 2 and \t\ < 1 GeV 2 , where xjp is the fractional energy loss of the final state proton 
and t the invariant momentum transfer, t = (P — P') 2 . The presence of large and positive scale 
violations up to the largest parton fractional momentum, /3, accessible to experiments reveals 
that DPDF's are gluon dominated distributions. Diffractive gluon distributions are however 
hardly determined by the diffractive structure functions alone and this fact has stimulated 
experimental collaborations to measure quantity that are directly sensitive to it El [9] and to 
include such data in global fits. 

The general method used to extract DPDF's from available data heavily relies on a number of 
assumptions motivated by Regge phenomenology. The purpose of this short note is to present 
a new method which does not require any Regge assumptions and is much closer in spirit to the 
factorization theorem. 



2. Data set and observable 



For this analysis we will use HI [TO] data in which DIS diffractive events, ep — > eXY , are selected 
by requiring a large rapidity gap between the hadronic system X and the low mass dissociative 
system Y (which may consist in the scattered p only). The kinematics of diffractive events is 
specified by the following variables: 

Q z + p sxb 

where Q 2 is the photon virtuality, M x is the invariant mass of the hadronic final state X, xb 
is the usual x-Bjorken DIS variable and -y/s = 318 GeV is the HERA center of mass energy. 
The measurament is integrated over the hadronic final state Y mass region M Y < 1.6 GeV and 
the invariant momentum transfer |t| < lGeV 2 , which then define the extracted DPDF's |llj . 
The kinematical coverage is wide ranging from 3.5 < Q 2 < 1600GeV 2 , 3 • 10 -4 < xjp < 3 • 1CP 2 
and 10 -3 < P < 0.8. We notice however that the {/3,Q 2 } coverage changes with xjp due to 
kinematic constraints and this fact will reflect on DPDF's extraction. Data are presented as a 
three-fold reduced e + p cross section which depends on the diffractive structure functions 
and i*£ 3 ■ In the one-photon exchange approximation, it reads: 

a^^,Q 2 ,xjp) = F 2 D ^^,Q 2 ,xjp) - 1 + ^_ y)2 F^ 3) ((3,Q 2 ,xjp). (2) 



3. The new method 



The widely used approach \TU[ [T2| [T3] to extract DPDF's is to assume proton vertex 
factorization, i.e. that DPDF's can be factorized into a flux factor depending only on xp 
and t and a term depending only on /3 and Q 2 : 

ff 09, Q 2 , xjp, t) = fjp /P (xjp, t) Ff(/3, Q 2 ) + fn/pixjp, t) Q 2 ) + ... (3) 

Each term in the expansion, according to Regge theory, is supposed to give a dominant 
contribution in a given range of xjp, the pomeron (IP) at low xjp and the reggeon (M) at 
higher value of xjp. The flux factor fjp/p (fip/p) can be interpreted as the probability that a 
pomeron (reggeon) with a given value of xjp and t couples to the proton. The fluxes xjp and t 
dependences used in the fits are motivated by Regge theory. This approach requires that parton 
distributions of the pomeron, Ff, and of the reggeon, Ff 1 , must be simultaneously extracted 
from data. This procedure introduces a large number of parameters in the fit and it is potentially 
biased by the choices of the flux factors. A common choice used in phenomenological application 
is to fix F^ to be equal to pion parton distribution functions. Although such an approach has 
been proven to be supported by phenomenological analyses within the precision obtained with 
HERA-I data, it is not routed in perturbative QCD and may show up to be not satisfactory 
with the expected precision increase of HERA-II data and Hl+ZEUS complete combination. 
The alternative method we propose is instead inspired by the factorization theorem [1] for 
diffractive DIS itself. The latter states that factorization holds at fixed values of xjp and t so 
that the parton content described by Ff is uniquely fixed by the kinematics of the outgoing 
proton and it is in principle different for different values of xjp. In practice this idea is realized 
performing a series of separate pQCD fits at fixed values of xjp with a common initial condition 
controlled by a set of parameters {pi}- This procedure guide us to infere the approximate 
dependence of parameters {p{\ on xjp allowing the construction of initial condition in the {f3, xjp} 
space to be used in a global fit, without any further model dependent assumption. 



Table 1. Fit results at fixed xpom. Fits are performed taking into account only statistical 
errors only. Only data for which M 2 > 4 GeV 2 and Q 2 > 8.5 GeV 2 are included in the fits. 



Xjp 


X 2 1 d.o.f 


0.001 


0.921 


0.003 


0.875 


0.01 


0.882 


0.03 


0.472 



4. Fit procedure and results 

In each xjp-bin for which data are presented by the experimental collaboration we perform a 
separate pQCD fits. For the singlet and gluon distributions at the arbitrary scale Q\ we choose: 

/3S(/3,Q 2 ) = A g p B *(l-/3f«e-^, 

(3g(f3,Q 2 ) = A g e-°^, (4) 

so that there are four free parameters. We make the common assumption that all lights quark 
distributions are equal to each other and the exponential dumping exponential factor allows 
more freedom in the variation of the parameters C q at large j3. The functional form in eq. @ 
is identical to the one used by HI collaboration in Ref. [10] . Since diffractive DIS data hardly 
discriminate [10J between different behaviour of gluon distributions at large (3 we choose the 
simpler one in which the gluon is a costant at Qq. Such distributions are then evolved with 
the QCDNUM17 [2] program within a fixed flavour number scheme to next-to-leading order 
accuracy. Heavy flavours contributions are taken into account in the general massive scheme. 
The convolution engine of QCDNUM17 is used to obtain F^ S ' and fF^ structure functions at 
next-to-leading order which are then minimized versus data. In order to avoid the resonance 
region, a cut on the invariant mass of the hadronic system X is applied, M 2 > 4 GeV 2 . Single 
xjp-fit results are sensitive to the choice of the mininum value Q 2 of data to be included in 
the fits. The inclusion of data with Q 2 < 8.5 GeV 2 in general worsen the x 2 an d induce large 
fluctuation in the gluon distribution. This instability has been already noticed in Ref. |10| and 
avoided by including in the fit only data for which Q 2 > 8.5 GeV 2 . We will adopt here the same 
strategy. 

An essential condition for the procedure to work is that good quality fits are all obtained with 
the common initial condition, eq. Q. From fits results presented in Table [1] we conclude that 
the initial condition provided by eq. is general enough to describe data in all xjp-bins. The 
dependence of the parameters on xjp is reported in Figures from 1 to 4. Red dots are the results 
from pQCD fits at fixed xjp. The singlet normalization A q behaves as an inverse power of xjp. 
In order to improve the description at higher xjp, however, an additional term is also included: 

A q {x P ) = A qfi (xip) a ^ (1 - x F ) A ^ . (5) 

The gluon normalization is compatible with a single inverse power behaviour of the type: 



A g { X] p) = A gfi (xjp) a °' 



(6) 
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Figure 1. A q as a function of xjp. 



Figure 2. A g as a function of xjp. 




0.005 0.01 0.015 0.02 0.025 0.O3 0.035 




11.005 0.01 



0.015 0.02 0.025 0.O3 0.035 



Figure 3. B q as a function of xjp. 



Figure 4. C q as a function of xjp. 



The coefficients B q and C q which control the /3-shape of the singlet distribution are well described 
by: 



B q (xip) 



Bqfl + -Bg,l X]p , 
Cqfi + Cq,l X]P ■ 



(7) 
(8) 



In order to facilitate the comparison, curves resulting from fits with functional form specified 
in eqs. (6-9) are superimposed to points in Figures 1-4. Since the approximate behaviour of 
the parameters against xjp is now known, we can perform a a^-bin combined fit, using the 
generalized initial condition 



/3S(/3,Q2 iX7p ) 

P g(/3,Qo,xip) 



A q { XlP ) (3 B "^ (1-/3) 

0.01 

A g (x]p) e , 



(9) 
(10) 



with the parameters depending on xjp as specified in eqs. (6-9). The combined fit has nine free 
parameters. Following the procedure described in Ref. [15], to each systematic errors quoted 
in the experimental paper is assigned a free systematic parameters which is then minimized in 
the fit along with parameters. As for the single-xjp fits, only data for which > 4GeV 2 and 
Q 2 > 8.5 GeV 2 are included in the fit. The latter has a appreciable sensitivity on the scale Qq 
due to the relative stiffness of the initial condition in eq. . The choice of Qq is then optimized 
performing a scan which gives the best x 2 value for Qq = 2.3 GeV 2 . The partial results for the 
X 2 function in the various xp bin are reported in Table [2j The best fit returns a \ 2 = 166 for 182 
degrees of freedom which is of comparable quality as the one presented in Ref. [10] . The initial 
condition, eq. ©, allows the singlet and gluon normalization, A q and A g respectively, to have a 
different power behaviour. It is therefore interesting to notice that if the condition A qt i = A gt \ 
is imposed, this result in a global x 2 = 171 for 183 de gree of freedom. If one further neglects 



Table 2. Combined xp fit result. 



X]p 


x 2 


fitted 


total 






points 


points 


0.0003 


3.1 


3 


12 


0.0010 


26.2 


19 


40 


0.0030 


38.9 


39 


59 


0.0100 


58.6 


61 


80 


0.0300 


36.4 


68 


85 



the xjp-dependence of B q and C g by setting B q ^\ = C ?il = the x 2 increases to 188 units for 
185 degree of freedom. This is an a posteriori confirmation that not only diffractive parton 
distributions change their magnitude versus x jp but also that a modulation in their /3-shape (for 
the singlet, in this case) is necessary to better fit the data. 



5. Conclusions 



We have outlined a new method to extract diffractive parton distributions inspired by the 
factorization theorem for diffractive DIS. From a series of pQCD fits at fixed sp we infere the 
dependence of parameters on such a variable and this allows us to construct a generalized initial 
condition without assuming neither proton vertex factorization nor the existence of a series of 
Regge trajectories. Although the quality of the resulting fit gives a \ 2 I d.o.f close to unity, as 
the Regge based pQCD fit [10J, the new procedure treats the non-perturbative xjp-dependence of 
the cross-sections in a controlled and less model dependent way. This feature will allow further 
factorization test to be performed with the improved precision of published [16] and forthcoming 
data OH IB]. 
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